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Traditionally, the quantum Brownian motion is described by Fokker-Planck or diffusion equations 
in terms of quasi-probability distribution functions, e.g., Wigner functions. These often become 
singular or negative in the full quantum regime. In this paper a simple approach to non-Markovian 
theory of quantum Brownian motion using true probability distribution functions is presented. Based 
on an initial coherent state representation of the bath oscillators and an equilibrium canonical 
distribution of the quantum mechanical mean values of their co-ordinates and momenta we derive 
a generalized quantum Langevin equation in c-numbers and show that the latter is amenable to a 
theoretical analysis in terms of the classical theory of non-Markovian dynamics. The corresponding 
Fokker-Planck, diffusion and the Smoluchowski equations are the exact quantum analogues of their 
classical counterparts. The present work is independent of path integral techniques. The theory as 
developed here is a natural extension of its classical version and is valid for arbitrary temperature 
and friction (Smoluchowski equation being considered in the overdamped limit). 

PACS number(s) : 05.40.-a, 05.30. Ch, 02.50.-r 
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O , I. INTRODUCTION 

A model quantum system coupled to its environment forms the standard paradigm of quantum Brownian motion. 
The initiation of early development of this stochastic process took place around the middle of this century A 
, major impetus was the discovery of laser in sixties followed by significant advancement in the field of quantum optics 
and laser physics in seventies where the extensive applications of nonequilibrium quantum statistical methods were 
made. Various nonlinear optical processes/phenomena were described with the help of operator Langevin equations, 
density operator methods and the associated quasi-classical distribution functions of Wigner, Glauber, Sudarshan and 
others centering around the quantum Markov processes jj]-^). Subsequent to this early development the quantum 
theory of Brownian motion again emerged as a subject of immense interest in early eighties when the problem of 
macroscopic quantum tunneling was addressed by Leggett and others ||7|— [1 1|| and almost simultaneously quantum 
Kramers' problem attracted serious attention of a number of workers p^-|l5|. The method which received major 
appreciation in eighties and nineties in the wide community of physicists and chemists in these studies is the real time 
functional integral |l6| , |l7[ |. This method has been shown to be an effective tool for treatment of quantum transition 
state [pH , dissipative quantum coherence effects 19 as well as incoherent quantum tunneling processes ]b^[i4] , pc| 



and many related problems pi, , . 

Inspite of this phenomenal success it may, however, noted that compared to classical theory quantum theory of 
Brownian motion based on functional integrals rests on a fundamentally different footing. While the classical theory 
is based on the differential equations for evolution of true probability density functions of the particle executing 
Brownian motion, the path integral methods rely on non-canonical quantization procedure and the evaluation of 
quantum partition function of the particle interacting with the heat bath and one is, in general, led to the time 
evolution equations of quasi-probability distribution functions such as Wigner functions [fb5| , p2| |2(f| . The question is 
whether there is any natural extension of classical method to quantum domain in terms of true probability distribution 
functions. It is therefore worthwhile to seek for a natural extension of the classical theory of Brownian motion 
to quantum domain in the non-Markovian regime for arbitrary friction and temperature within the framework of a 
well-behaved true probabilistic description. 
Our aim in this paper is thus twofold : 

(i) to enquire whether there exists a quantum generalized Langevin equation (QGLE) in c-numbers whose noise 
correlation satisfies the quantum fluctuation-dissipation relation (FDR) but which (QGLE) at the same time is a 
natural analogue of its classical counterpart. 

(ii) to formulate the exact quantum Fokker-Planck and diffusion equations which are valid for arbitrary temperature 
and friction. We also intend to look for the overdamped limit to obtain the exact quantum analogue of classical 
Smoluchowski equation. 
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Before proceeding further it is important to stress the motivation for the present scheme: 

(1) As we have already pointed out that the traditional theories of quantum Brownian motion in optics and 
condensed matter physics ^] are based on quasi-probability functions. Apart from their usual shortcomings that they 
may become negative or singular |27j| in the full quantum regime when the potential is nonlinear, the quasi-probability 
functions are, in general, not valid for non-Markovian processes with arbitrary noise correlation. While in majority 
of the quantum optical situations Markovian description is sufficient, non-Markovian effects of noise correlation are 
strongly felt in the problem of quantum dissipation in condensed matter and chemical physics at low temperature. 
To include these effects even in the case of a free particle [ see for example, Ref. jll] ] one has to use a suitable 
cut-off frequency of the heat bath to avoid intrinsic low frequency divergence. Clearly this poses serious difficulties for 
studying transient behavior for arbitrary noise correlation and temperature. In what follows we show that the present 
treatment is free from such difficulties. 

(2) Our second motivation is to understand quantum- classical correspondence in the problem of Brownian motion 
in a transparent way. To this end we note that in the classical theory the Fokker-Planck equation with nonlinear 
potential contains derivatives of probability distribution functions upto second order. The equations in terms of 
Wigncr functions on the other hand involve higher (than two) order derivatives of distribution functions in the 
corresponding quantum formulation ps|] . The higher derivative terms contain powers of h and derivatives of potential 
signifying purely quantum diffusion in which quantum corrections and nonlinearity of the potential get entangled in 
the description of the system. Because of the occurence of higher derivatives the positivity of the distribution function 
is never ensured and the equation cannot be treated as a quantum analogue of classical Fokker-Planck equation. Any 
attempt to reduce the order of the derivatives to two amounts to a semiclassical approximation. Again there exists 
no systematic procedure for this reduction. Keeping in view of these problems we intend to derive exact quantum 
analogues [ Eqs.(42), (47) and (50) ] of classical Fokker-Planck, diffusion and Smoluchowski equations, respectively, 
in terms of true probability distribution function where the equations contain derivatives of distribution functions 
upto second order only for which the diffusion coefficients are positive definite. Since the equations are classical 
looking in form but quantum mechanical in their content one can read the quantum drift and diffusion coefficients 
and also construct the quantum corrections due to the nonlinearity of the system systematically order by order in a 
straightforward way so that quantum-classical correspondance can be checked simply by taking limit h — * both in 
Markovian and in non-Markovian description. We mention in passing that in contrast to a recent treatment [^9| of 
large friction limit in a similar context, the quantum Smoluchowski equation as discussed here retains its validity in 
the full quantum regime as T — ► 0. 

(3) Since over the last two decades classical non-Markovian theories 30 3l]] and numerical methods of generating 
classical noise processes have made a significant progress |j32|-|34| , the mapping of quantum theory of Brownian motion 
into a classical form, as achieved here, suggests that the classical treatment can be extended to quantum domain 
without much difficulty. Since the present scheme describes the generation of quantum noise [ Eqs. (10) and (If) ] 
as classical numbers which follow quantum fluctuation-dissipation relation it is easy to comprehend that the classical 
numerical techniques of generation of noise and solving stochastic Langevin equation can be utilized in the 
present case in a straightforward way to solve quantum Langevin equation [|35f . The procedure is therefore much easy 
to implement compared to other methods like path integral Monte Carlo techniques |36|. 

In what follows we consider the standard system-reservoir model and make use of the coherent state representation 
of the bath oscillators to derive a GLE for quantum mechanical mean value of position of a particle in contact with 
a thermal bath whose quantum mechanical properties can be defined in terms of a classical-looking noise term and a 
canonical distribution of initial quantum mechanical mean values of the co-ordinates and momenta of the bath. This 
simple approach allows us to show that although the equation is essentially quantum mechanical it is amenable to a 
theoretical analysis in terms of the classical theory of non-Markovian dynamics |3(],|3l[] . 

The rest of the paper is organized as follows : The system reservoir model, the associated QGLE and the canonical 
distribution for the bath oscillators have been introduced in Sec. II. This is followed by a general analysis of QGLE 
in Sec. Ill and an illustration with an exponential memory kernel in Sec. IV to calculate the variances required for 
setting up a quantum Fokker-Planck equation and a quantum diffusion equation in sections V and VI, repectively. 
Section VII is devoted to quantum overdamped limit and Smoluchowski equation. The paper is summarized and 
concluded in Sec. VIII. 



II. THE QUANTUM GENERALIZED LANGEVIN EQUATION (QGLE) IN C-NUMBERS 

We consider a particle in a medium. The latter is modeled as a set of harmonic oscillators with frequency {tOi}. 
Evolution of such a quantum open system has been studied over the last several decades under a variety of reasonable 
assumptions. Specifically our interest here is to develop an exact description of quantum Brownian motion within the 
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perview of this model described by the following Hamiltonian 



H = 



1 



(1) 



Here x and p are co-ordinate and momentum operators of the particle and the set {qj,Pj} is the set of co-ordinate and 
momentum operators for the reservoir oscillators coupled linearly to the system through their coupling coefficients 
Kj. The potential V(x) is due to the external force field for the Brownian particle. The co-ordinate and momentum 
operators follow the usual commutation relation [x,p] = ih and [qj,pj] — ihdij. Note that in writing down the 
Hamiltonian no rotating wave approximation has been used. 

Eliminating the reservoir degrees of freedom in the usual way p|,|3^ -flofi we obtain the operator Langevin equation 
for the particle, 



x(t) + / d£ j(t - t') x(t') + V'{x) = F{t) 



(2) 



where the noise operator F(t) and the memory kernel j(t) are given by 



F(t) = [ {?i(0) - i(0)} cosuijt + pj(0) k 1 / 2 sin 



LVjt 



(3) 



and 



7(i) = Kj cosujjt 



(4) 



with m 



( masses have been assumed to be unity 



The Eq.(|^) is an exact quantized operator Langevin equation which is now a standard textbook material 
and for which the noise properties of F(t) can be defined using a suitable initial canonical distribution of the bath 
co-ordinates and momenta. Our aim here is to replace it by an equivalent QGLE in c-numbcrs. Again this is not a 
new problem so long as one is restricted to standard quasi-probabilistic methods using, for example, Wigner functions 
fl5| , p2 - 26 1 . To address the problem of quantum non-Markovian dynamics in terms of a true probabilistic description 
we, however, follow a different procedure. We first carry out the quantum mechanical average of Eq.(Q) 



(*(«)) 



dt><y(t-t>) (x(t')) + {V'(x))=(F(t)) 



(5) 



where the average (. . . 
at t = 0, |0){|ai)|a 2 ) 



is taken over the initial product separable quantum states of the particle and the bath oscillators 
. . . |ajv)}' Here \<p) denotes any arbitrary initial state of the particle and jo^) corresponds to 



the initial coherent state of the i-th bath oscillator. \on) is given by \on) = exp(— |a<j| /2) J2^ i= o( a T / v™^)\ n i) > 
cti being expressed in terms of the mean values of the co-ordinate and momentum of the i-th oscillator, (<7i(0)) = 
(Vh/2coi)(ai + a*) and (pi(0)) = iy/hidi/2(a* — cti), respectively. It is important to note that (F(t)) of Eq.(||) is a 
classical-like noise term which, in general, is a non-zero number because of the quantum mechanical averaging over the 
co-ordinate and momentum operators of the bath oscillators with respect to the initial coherent states and arbitrary 
initial state of the particle and is given by 



(F(t)) = [ {fe(0)) - (x(0))} Kj cosujt+(pj(0)) k) 
j 

It is convenient to rewrite the c-number equation as follows ; 



/2 



sincjjt 



(6) 



(x(t)) + / dt' 7 (t - t') (x(t')) + (V'(x)) = F(t) 
Jo 



(7) 



where we let the quantum mechanical mean value (F(t)} = F(t). We now turn to the second averaging. To realize F(t) 
as an effective c-number noise we now assume that the momenta (j)j(0)) and the shifted co-ordinates {(<2j (0)) — (x(0))} 
of the bath oscillators are distributed according to a canonical distribution of Gaussian forms as 
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Vj = TV exp • 



[fe(o)) 2 + ^{(*(o))-(i(o))} 2 



2hUj (nj + i) j 

so that for any quantum mechanical mean value Oj((pj(0)) , {(qj(0)) — (x(0))}) the statistical average (. . is 



(8) 



(Oj)s = J O,(fe(0)),{fe(0))-(i(0))}) 

x ^-((^(0)>,{(^(0)> - (*(0)») d(p 3 (0)) d{(ft(0)) - (x(0))} 



(9) 



Here fij indicates the average thermal photon number of the j-th oscillator at temperature T and fij = 
l/[exp (hujj /ksT) — 1] and Af is the normalization constant. 
The distribution (||) and the definition of statistical average 



imply that F(t) must satisfy 







and 



(F(t)F(t')) s 



2 ^ 



(10) 



(11) 



That is, the c-number noise F(t) is such that it is zero centered and satisfies the standard quantum fluctuation- 
dissipation relation (FDR) as known in the literature m terms of quantum statistical average of the noise operators. 
To proceed further we now add the force term V'{(x)) on both sides of Eq.(^|) and rearrange it to obtain formally 



X(t) + / dt' 7 (f - t') X(t') + V'(X) = F(t) + Q(X, t) 
Jo 

where we let (x(t)) = X(t) for simple notational convenience and 

Q(X,t) = V'({x))-{V'(x)) 



(12) 



(13) 



represents the quantum mechanical dispersion of the force operator V'(x) due to the system degree of freedom. Since 
Q(t) is a quantum fluctuation term Eq.(12) offers a simple interpretation. This implies that the classical looking 
QGLE is governed by a c-number quantum noise F(t) which originates from the quantum mechanical heat bath 
characterized by the properties ( p~0| ) and ( |ll| ) and a quantum fluctuation term Q(t) due to the quantum nature of the 
system characteristic of the nonlinearity of the potential. In Sec. VII we give a recipe for calculation of Q(t). 

Summarizing the above discussions we point out that it is possible to formulate a QGLE (12) of the quantum 
mechanical mean value of position of a particle in a medium, provided the classical-like noise term F(t) satisfies (|l(i| ) 
and (|ll|) where the ensemble average has to be carried out with the distribution (||) . It is thus apparent that to realize 
F(t) as a noise term we have split up the standard quantum statistical averaging procedure into a quantum mechanical 
mean (. . .) by explicitly using an initial coherent state representation of the bath oscillators and then a statistical 
average (■ ■ -)s of the quantum mechanical mean values. Two pertinent points are to be noted : First, it may be 
easily verified that the distribution of quantum mechanical mean values of the bath oscillators (^) reduces to classical 
Maxwell-Boltzmann distribution in the thermal limit, hujj <C ksT. Second, the vacuum term in the distribution (j^) 
prevents the distribution of quantum mechanical mean values from being singular at T = ; or in other words the 
width of distribution remains finite even at absolute zero, which is a simple consequence of uncertainty principle. 



III. GENERAL ANALYSIS: DAMPED FREE PARTICLE 

It is now convenient to rewrite QGLE (|l^) of quantum mechanical mean value of position of a particle in the 
absence of any external force field in the form 

X(t) + [ j(t- t') X(t') dt' = F(t) (14) 



7(i) is the dissipative memory kernel as given by Eq.(|) and F(t) is the zero centered stationary noise, i.e., 

(F(t)) s = and (F(t)F(t')) s = C(\t - t'\) = C(r) (15) 
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where C(t) is the correlation function which in the equilibrium state is connected to the memory kernel j(t) through 
FDR of the form M 



1 f°° ( ?ioj \ 

C(t — t') = — / duj k(oj) q(uj) Tiuj ( coth I cosw(i — i') 

2 Jo " V 2fc B T/ 



(16) 



Eq.(|16[) is the continuum version of Eq.(|ll|). q(oj) denotes the density of modes of the bath oscillators. Here it is 
important to note that Eq.(|l6|) is the generalized FDR valid at any arbitrary temperature T. ^{t—t') is the continuum 
version of Eq. (Q) and is given by 



j(t-t')= / duo h(uj)q(uo) cos uj(t-t') . (17) 
Jo 

In the high temperature limit, i.e., for hto <C fc^T we arrive at the wellknown classical FDR of the second kind Jll]] 

C(t - t') = k B T 7 (f - t') . (18) 
The general solution of Eq.(^i|) is given by 

X(t) = (X(t)) s + I H(t-T)F(r)dT (19) 
Jo 

where 

(X(t)) s = X Q + V H(t) (20) 

with Xq = X(0) and Vq — X(0) being the initial quantum mechanical mean values of position and velocity of the 
particle, respectively. H (t) is the inverse form of the Laplace transform 

- 2 * , v (21) 

with 

/•OC 

j(s) = / j^e-^dt (22) 
Jo 

is the Laplace transform of dissipative memory kernel 7(t). The time derivative of Eq.(|l9|) gives 

V(t) = (V(t)) s + I h(t- t) F(t) dr (23) 



where 

(V(t)) s = V h(t) (24) 

and 

m - « ■ <*) 

Hence 

Ms) = — • ( 26 ) 

Before proceeding further it is important to recall the physical significance of the two function H(t) and h(t). It 
has already been assumed that the initial quantum mechanical velocity Vb is independent of the random force F(t), 

(V F(t)) s = . (27) 
Thus multiplying Eqs. © and (||) by V and using relation ^27^) we obtain, 
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(V V(t))s/(V Q 2 ) s = h(t) 



(28) 



(V (X(t)-Xo))s/(V 2 ] 



H(t) 



(29) 



Hence H{t) and h(t) are the two relaxation functions ; h(t) measures how the quantum mechanical mean velocity 
forgets its initial value and H (t) measures how the quantum mechanical mean displacement forgets the initial velocity. 
As a result quantum mechanical mean velocity of the particle relaxes to a stationary state with zero statistical average 
of the quantum mechanical mean velocity. 

Now using the symmetry property of the correlation function 

(F(t)F(t')) s = C(t-t')=C(t'-t) 

and using the solution for X(t) and V(t) we obtain the following expressions of the variances, 



'xx 



>vv 



(t) = ([X(t)-(X(t)), 



(t) 



2 f H{h) dh f 1 H(t 2 ) C(h - t 2 )dt 2 
Jo Jo 

(lV(t)-(V(t)) s f) s 



h(ti) dh / h(t 2 ) C(h 



t 2 )dt 2 and 

([X(t) (X(t)) s ][V(t) (V(t)) s ])s = \*xx{t) 



H(h) dh / h(t 2 ) C(h - t 2 )dt 2 



(30a) 



(30b) 



(30c) 



The above three expressions are valid for arbitrary temperature and friction and include quantum effects. However 
in the high temperature classical limit ( i.e., Tiuj <C ksT ) one can derive simplified versions of the variances 



a\ x (t) = k B T 2 f H(t')dt' - H 2 (t) 
I Jo 

iy V (t) = k B T[l-h 2 {t)] and 
aj cv {t)=k B TH{t)[l-h(t)] . 



(31a) 

(31b) 
(31c) 



Before closing this section we emphasize a pertinent point at this stage. The (30a)-( p0q ) are the expressions for 
statistical variances of the quantum mechanical mean values X and V. These are not to be confused with the standard 
quantum mechanical variances which are connected through uncertainty relations. 



IV. A SPECIFIC EXAMPLE : EXPONENTIALLY CORRELATED MEMORY KERNEL 



The very structure of j(t) given in Eq.(|17]) suggests that it is quite general and a further calculation requires a 
prior knowledge of the density of modes g(u>) of the bath oscillators. As an specific case we consider in the continuum 
limit. 



k(u>)q(u>) = — 



7o 



so that 7(i) takes the wellknown form, 



lit) 



n 1 + ui 2 t^ 



Tr 



(32) 



(33) 



where 70 is the damping constant and r c refers to correlation time of the noise. Once we get an explicit expression of 
j(t) in closed form and its Laplace transform, it is possible to make use of Eq. (pl|) to calculate the relaxation function 
H(t), which for the present case is given by 



H(t) = — [l - AT*/ 2t ° sin(At + a) 
7o L 



(34) 



G 



where 



A V T r . 4t? 



1/2 



and a — tan 



2Ar r 



1 - 2-f T c 



(35) 



Now making use of the expressions for H(t) and the correlation function C(t) in Eqs.(30a)-(30c) we calculate 
explicitly after a long but straightforward algebra the time dependent expressions of the variances of the quantum 
mechanical mean value of position and momentum of the particle, 



'xx 



* = — / 



coth 



2k B Tj 



Tx(y>, t) dus 



(36) 



and 

<4v(*) = ^xxW (38) 

In Appendix-A we provide the explicit structures of Tx{oJ,t) and fviujt). 

To examine the consistency of our calculation we check long time behaviour of the classical high temperature Ohmic 
limit of the variances u xx (t) and <Jy V (t) . In this limit we have 

2 . , 4k B T f°° 1 . . 

<Txx{t) = —T~ I duJ -,,,22 Sx(U,t) 







1 + L0 2 T r 2 



Only the first term of Tx{oJ,t) gives the long time behaviour of <J xx (t) in the Markovian limit, contribution of the 
rest of the terms being zero. Taking this leading order contribution we have 



Ak B T r°° ^ 1 1 

7T7o J 1 + W 2 T C 2 OJ 

8k B T ( 1 



a xx(t) = -z— / du - ; |2 ^ 2 — (1 - cosujt) 

sin 2 i^i 
1 — 



7T7o V 1 + w 2 r 2 
which gives 



u=0/ JO 



4xW = - — t fori^cx) . (39) 



To 

Similarly one can show that for classical high temperature Markovian limit 

crlr V (t)^k B T fort^cx) . (40) 

Since we are unable to evaluate analytically further the explicit time dependent structures of the variances in the 
general case, we take resort to numerical integration of Eqs.(^) and (p7|) . In Figs.(l) and (2) we show the short 
time and long time behaviour of the variances o xx (t) as functions of time for different values of temperature but for 
a fixed value of correlation time, r c . It is apparent that while the short time dynamics has a simple t 2 behaviour, 
asymptotic dependence is linear in t with a clear cross-over around some intermediate time. Fig. (3) exhibits the 
asymptotic constancy of Oy V (t) as a function of time for different temperatures. The effect of correlation time r c on 
the variance (?xx(t) nas been examined in Fig. (4) for a fixed high temperature k B T = 10.0. It is interesting to note 
that the cross-over region gets longer for larger correlation time. 

Figs. (5) and (6) illustrate the zero temperature situation. In this regime non-Markovian effects are strong which 
is evident from vacuum fluctuations growing in time in an oscillatory fashion at early stages for different values of 
correlation time as shown in Fig. (5). In Fig. (6) we show how the initial growth of variance <Jy V (t) finally settles down 
to a constant non-thermal energy value. 
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V. THE GENERALIZED QUANTUM FOKKER-PLANCK EQUATION 



We now return to our general analysis as carried out in Sec. III. To write down the Fokker- Planck description for 
the evolution of probability density function of quantum mechanical mean values of co-ordinate and momentum of 
the particle it is necessary to consider the statistical distribution of noise which we assume here to be Gaussian. For 
Gaussian noise processes we define the joint characteristic function in terms of the standard mean values and variances 
as follows : 



p {(i>,P,t) = exp 



ilM{X(t)) s + ip(V(t)) s - - {o 2 xx {t)p? + 



(41) 



Using the standard procedure |3(]j3l|] we write down below the Fokker-Planck equation (FPE) obeyed by the joint 
probability density function P(X, V, t) which is the inverse Fourier transform of the characteristic function : 

di + v dx) p[x ' v ' t] = m dv VP[X ' v ' t] + ^ dv^ p(x ' v > l) 

^dxW^^ (42) 

where 

at) = -h(t)/h(t) , (43a) 

¥>(*) = £(t)vvv{t) + ~dvv(t) and ( 43b ) 

m = -°vv(t) + Z(t>xv(t) + °xv(t) ■ (43c) 

The above FPE is the exact quantum mechanical version of the classical non-Markovian FPE and is valid at any 
arbitrary temperature and friction. 

The decisive advantage of the present approach is again noteworthy. We have mapped the operator generalized 
Langevin equation into a generalized Langevin equation in c- numbers (|l4) and its equivalent Fokker-Planck equa- 
tion ([42]). The present approach bypasses the earlier methods of quasi-probabilistic distribution functions employed 
widely in quantum optics over the decades in a number of ways. First, unlike the quasi-probabilistic distri- 

bution functions, the probability distribution function P(X, V, t) is valid for non-Markov processes. Second, while 
the corresponding characteristic functions for quasi-probabilistic distribution functions are operators, we make use 
of characteristic functions which are numbers. Third, as pointed out earlier the quasi-distribution functions often 
become negative or singular in the strong quantum domain and pose serious problems. The present approach is free 
from such shortcomings since the probability density function, P{X, V, t) behaves here as a true probability function 
rather than a quasi-probability function. 



VI. GENERALIZED QUANTUM DIFFUSION EQUATION 

In their landmark paper on classical Brownian motion Ornstein and Uhlenbeck |12| solved the classical Markovian 
FPE to find P(X, V, t) and then in a bid to obtain Einstein's diffusion equation tried to evaluate p(X, t), the probability 
density function in configuration space by integrating over V . It was shown that it is difficult, if not impossible to 
obtain a differential equation for P(X,Vo,t) from the classical Markovian FPE which for t ^> I/70 would become a 
diffusion equation. However, for the classical non-Markovian case Mazo |[o| in late seventies addressed this problem 
by considering an initial Maxwellian distribution $(Vb) of the initial velocity Vb and then derived the exact differential 
equation satisfied by p{X, t) where 

p(X,t) = J P(X,V ,t) $(Vb) dV . 

The resulting equation thus reduces to the diffusion equation for t 3> 1/7. We follow Mazo's procedure to derive an 
exact quantum mechanical version of the classical non-Markovian case, a differential equation which for t > 1/7 goes 
over into a quantum diffusion equation. To this end we proceed as follows; from Eq. (|4l|) for p — case we have 



8 



p(ji,t) = J p(/*,i) nVo)dV 

= exp (-~M 2 <4 x (t)) exp(^Xo) | exp[i/iV ff(t)] *W>) ■ (44) 
Here we take the initial Gaussian distribution of the quantum mechanical mean values of the velocity of the particle, 

*(Vo)=(^r-) 12 exp(S) (45) 



V27rA o y V 2A 

where 

A = ¥>(oo)/£(oo) . (46) 
It is not difficult to note that the above choice is dictated by the stationary solution of the QFPE (|42|) 



satisfies (|42|) at equilibrium. The explicit time-dependent expressions for ip(t) and £(t) have been given in (43a) and 



(43t ). Inserting Eq.(ktq) in (nM and then performing the inverse Fourier transform after integration over Vq we arrive 



at the following equation after little a algebra, 

dt q[ ) ax* ■ [ t) 

This is the quantum analogue of Einstein's diffusion equation where the explicit structure of the time-dependent 
quantum diffusion coefficient, D q (t) is given by, 

D q (t)=a xv (t)+A H(t)h(t) . (48) 



The required variances, the relaxation functions and other related quantities in Eq.(|48|) are given in (30c), (|2^), ( |2l| ) 
and Jig). We now discuss the limiting cases. For classical Markovian limit the variance <J xv (t) gives fcsT/7o for 
t ^> I/70 and the second term in D q (t) vanishes in the long time limit, so that we recover Einstein's diffusion coefficient 
in configuration space. In the low temperature, however, the quantum effects begin to dominate. It is interesting 
to note that based on Feynman- Vernon path integral technique |f6|,[l7j , Hakim and Ambegaokar jllj had considered 
explicit quantum corrections to classical diffusion to examine the differential behavior of high and low temperature 
dependence in the dynamics for Leggett-Caldeira initial conditions. The non-Markovian nature of the dynamics is 
taken into account by considering the frequency dependence of the bath with a suitable low frequency cut-off. The 
transient behavior in the quantum correction to classical diffusion is therefore only observable on the timescales longer 
than the inverse cut-off frequency. The present treatment being exact, equipped to deal with arbitrary noise correlation 
at all temperatures and free from divergences does not require any such cut-off. The quantum diffusion coefficient can 
be followed arbitrarily from transient to the asymptotic regions. To explore the associated non-Markovian nature of 
the dynamics in the present case it is necessary to go over to numerical evaluation of D q (t). In Fig. (7) [ compare with 
Fig. I of Ref. |Tl[] ] we plot the variation of quantum diffusion coefficient D q {t) for several values of temperatures as 
a function of time for the exponential memory kernel considered in our example in Sec. IV. It is apparent that while 
the short time behaviour is characterized by a sharp increase followed by a maximum, the diffusion coefficient settles 
down to a constant value in the asymptotic limit. The short time behaviour is dominated by the second term in fl4S| ) 
due to the relaxation functions H (t) and h(t) of which the latter vanishes in the long time limit. Again the first term 
in ( p8^ offers no contribution to diffusion coefficient from its classical part in the vacuum limit at T = 0. The solid 
curve in Fig. (7) thus shows the evolution of a non-thermal diffusion coefficient of pure quantum origin. 



VII. QUANTUM SMOLUCHOWSKI EQUATION 

We now consider the diffusion of a particle in an external potential V(X) as described by QGLE (|f2|). In the 
overdamped limit we drop the inertial term X(t) and the damping kernel ^(t — t') is reduced to 70 S(t — t') for 
vanishing r c in (|35|). 70 is the Markovian limit of dissipation. Eq.(|f2[) then assumes the following form 

X + -[V'(X)-Q(X,t)} = ^ . (49) 
7o 7o 

Expressing V'(X) — Q(X,t) as a derivative of an effective quantum potential V qua nt(X,t) with respect to X, the 
equivalent description in terms of true probability distribution function p(X, t) is given by 



9 



dp(X.i) 1 d r , , . . N n <9 2 n 

= To dx M™ nt (x ' t)p{x ' • ( 50a ) 

with 

V; uant (X,t)=V'(X)-Q(X,t) (50b) 
where Q(X, i) is defined in (13). Here D qo is the quantum diffusion coefficient in the overdamped limit which can be 



obtained with the help of the following definition fll 



^ pt + At rt + At ^ 

2D qo = — dtt dt 2 —(F(t 1 )F(t 2 )) s . (51) 

At J t Jt 7o 

Here the correlation function (F (t\)F (fe))g/7o of the c-number quantum noise is given by Eq. ([!(]) in the continuum 
limit. We then make use of Eq.(p2|) for vanishing t c in (JsiJ) to obtain after explicit integration 

1 

27o' 

where the frequency uj in ( ]52] ) refers to linearized frequency of the nonlinear system |l| . We now discuss the classical 
and vacuum limits of the quantum Smoluchowski equation (50a). It is easy to check that in the limit Ti£j <C ksT, D qo 



D qo = — %ui[2n(u>) + 1] (52) 



reduces to Einstein's classical diffusion coefficient fcsT/70. At the same time Q(X,t) vanishes so that Vq Uant (X,t) 
goes over to V'(X) and one recovers the usual classical Smoluchowski equation. In the opposite limit as T — > 0, 
however, both quantum noise due to nonlinearity of the system and vacuum fluctuation orginating from the heat 
bath make significant contribution. D qo in this limit assumes the form fiQ/2'yo- In this context we refer to a recent 
treatment on large friction limit in quantum dissipative dynamics p9| to point out that the latter theory does not 
retain its full validity as T — > since the quantum noise of the heat bath disapp ears in the vacuum limit. 

The second noteworthy feature about the quantum Smoluchowski equation ( |50a|) is that unlike Wigner function 
based equations it does not contain higher order (higher than second) derivatives of p(X, t). The positive defi- 
niteness of the probability distribution function is thus ensured. 

It is i mportant to emphasize at this juncture that so far as the ge neral formulation of the theory is concerned, 



Eq.( p0a| ) contains quantum corrections to all orders. In this sense Eq.(50a) is formally an exact quantum analogue 
of classical Smoluchowski equation. To make it more explicit we return to the quantum mechanics of the system in 
Heisenberg picture to write the operators x and p as 

x(i) = (i(i)) + Sx and p{t) = (p{t)} + dp . (53) 

(x(t)) and (p{t)) are the quantities signifying quantum mechanical averages and Sx and Sp are quantum corrections. 
By construction (Sx) and (Sp) are zero and they obey the commutation relation [<Jx,<5p] = ih. Using (|53| ) in (V'(x)) 
and a Taylor expansion around (x) it is possible to express Q(X,t) as [ see Eq.(13) ]: 

Q(X,t) = -Y^-V n+1 (X)(5x n (t)) (54a) 
n! 

n>2 



where V n (X) is the n-th derivative of the potential at X(= (x)). Eq.(54a) suggests a simple expression for an effective 
potential V quant (X,t) as 

V quant (X,t) = V(X) + ^Vn{X)(Sx n (t)) (54b) 

n>2 

where the classical potential V(X) gets modified by the quantum corrections to all orders. To solve quantum Smolu- 
chowski equation it is therefore necessary to calculate (Sx 2 (t)), (5x 3 (t)), etc. To the lowest order (x) and (Sx 2 ) follow 
a coupled set of equations as given below 

j t {x)=(p) (55a) 
±(p)=-V\(x))-\v'"{(x))(te 2 ) (55b) 

^-(Sx 2 ) = (SxSp + SpSx) (55c) 
dt 

^(SxSp + SpSx) = 2(Sf) - 2V"((x))(Sx 2 ) (55d) 



dt 

dt 



^-(Sf) = -V" ((x))(SxSp + SpSx) . (55e) 
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The above set of equations can be derived [[l3| from the Heisenberg's equation of motion. If one is interested in the 
local dynamics around a point (say, at the botto m or to p of the potential well), the set o f equations get decoupled and 



it is easy to obtain simple analytic solutions of (55a)-(55e) for (x) and (5x ) for (54a). The higher order estimates 



(e.g., fourth order) of the quantum corrections can be obtained from the solutions of the equations of successive 
higher order derived earlier by Sundaram and Milonni pq] or otherwise [E3|. Since the quantum corrections due to 
the system are calculated by different sets of equations for succesive orders, the measure of accuracy of truncation can 
be understood easily. It is, therefore, obvious that the present scheme provides a simple, systematic and quantitative 
estimate of the mean field and other decorrelation methods on the basis of quantum-classical correspondence. 

VIII. CONCLUSIONS 

The main purpose of this paper is to enquire whether a stochastic differential equation in c-numbers in the form of a 
generalized Langevin equation and its corresponding Fokkcr-Planck equation and diffusion equation and Smoluchowski 
equation in terms of true probability functions are viable for description of non-Markovian quantum Brownian motion. 
Based on an initial coherent state representation of bath oscillators and an equilibrium distribution of quantum 
mechanical mean values of their co-ordinates and momenta, which satisfy the essential properties of the associated 
noise of the bath degrees of freedom, we derive a QGLE for quantum mechanical mean value of the position of the 
particle. The main conclusions of this study are the following : 

(i) Our QGLE ([l4]) is amenable to analysis in terms of the methods developed earlier for the treatment of classical 
non-Markovian theory of Brownian motion. 

(ii) The generalized Langevin equation (|l2|), th e cor responding Fokker-Planck equation ( |42| ) and the diffusion 



equation ( |47| ) and also the Smoluchowski equation (50a) are the exact quantum analogues of their classical versions 
ptj| , ^l| . The probability distribution functions as employed here bear the true notion of statistical probability rather 
than that of quasi-probability. 

(iii) The theory of quantum Brownian motion developed here is valid for arbitrary noise correlation and temperature 
and is free from divergences. 

(iv) The realization of noise as a classical-looking entity which satisfies quantum fluctuation-dissipation relationship 
(O) allows ourselves to envisage quantum Brownian motion as a natural extension of its classical conterpart. The 
method is based on canonical quantization procedure and makes no reference to path integral formulations. 

We conclude by mentioning that the method discussed here is promising for simple differential equation based 
approaches [|l5| to quantum activated processes, tunneling problems as shown elsewhere fill , quantum ratchet jiF 
and in problems relating to the motion in periodic fields [pM>2] and allied issues. 
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APPENDIX A: THE EXPLICIT FORMS OF Fx{u,T) AND Fv{u,T) 
Txi^i t) consists of eleven terms which are given below ; 

T x (lo, t) = ( W , t) + J&> («, t) + (to, t) + Tf (w, t) + jf (w, t) 

+Tf (to, t) + (w, t) + (w, t) + (co, t) + (to, t) 

+J$ 1 \w,t) . (Al) 



The explicit structures of J 7 ^ (to, t) (i = 1, . . . , 11) are given by 



T K x> (w,t) = — (1 - cos ut) , (A2) 
7o^ 
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T^ x ' (co, t) = — [cos(a + u)t) — cos a] — [cos(a — cot) — cos a] 

7oto ' 70C0 

AA^ AA^ 

3 - [sin (a + Lot) — sin a] H — [sin(a — cot) — sin a] 



7 w 



7 co 



AA\ 



[ e -t/2r c { sin ( A ^ + a ) + 2Xt c cos(Ai + a)} 



27o 2 

{sin a + 2At c cos a}] 



4 4 W) = - 



l--^[ e -'/ 2 ^ {cos(At + a) - 2Ar c sin(At + a)} 
^7o 



{cos a — 2At c sin a}] 



-t/r c 



87o 2 

{sin 2a + 2Ar c cos 2a}] 



{sin 2(Xt + a) + 2Ar c cos 2(Xt + a)} 



! 1 



e~ t/Tc + ^ {2Ar c sin 2(Xt + a) - cos 2(Ai + a)} 

47qt c 



-t/r c 



470 T c 



(2At c sin 2a — cos 2a) 



A 



- — [4" ) {e~* /2Tc (2T c (A - co) sin[a + (A - w)t] - cos[a + (A - u)t]) 

7 CO 

— 2t c (A — co) sin a + cos a} 

-4" } {e"* /2Tc (2r c (A + w) sin[a + (A + to)t] - cos[a + (A + iv)t}) 
— 2t c (X + lo) sin a + cos a}] , 



2 



A 2 A. 



70 



:l [4 W ){ e -*/2rc( 2Tc ( A _ w ) sin ( A _ ^ t _ CQS ( A _ ^ + 1} 



-A^ ) {e" t/2T '=(2T C (A + co) sin[2a + (A + - cos[2a + (A + uj)t]) 
— (2r c (A + co) sin 2a — cos 2a)}] , 



A2A ^ ) [4" ) {e' t/2Tc (2T c (A + co) sin(A + w)i - cos(A + u)t) + 1} 
7o 

-A t f ) {e- t/2T "(2T c (X - to) sin[2a + (A - co)*] - cos[2a + (A - 
— (2t c (A — to) sin 2a — cos 2a)}] , 



A 2 A (^0 

.Fl 10) (co, i) = - 5 [^"^e-^'Csinpa + (A + u)t] + 2r c (A + co) cos[2a + (A + w)t]) 
7o 

— (sin 2a + 2t c (A + co) cos 2a)} 

+4 w) {e~* /2Tc (sin(A - w)i + 2r c (A - co) cos(A - w)t) - 2t c (A - co)}] 
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A2 a{^) 

^ U) (w,t) - 5-[4 w) {e~* /2Te (sin[2a + (A - + 2r c (A - w) cos[2a + (A - w)t]) 

7o 

— (sin 2a + 2t c (A — uj) cos 2a)} 

+A<f ) {e~ t/2Tc (sm(\ + uj)t + 2r c (A + uj) cos(A + uj)t) - 2t c (A + uj)}] 



where 



Similarly we have 



with 



A^ - T 



4 U) = 2r c 2 



1 



+ 



1 



1 + 4t 2 (A-cj) 2 l + 4r c 2 (A + w) 2 
A — cj A — uj 



A {uj) = — 
3 1 + 4t c 2 (A-w) 2 



1 + 4t 2 (A-cj) 2 1 + 4t 2 (A + cj) 2 



(W ) 2t 2 (A-u;) 
5 1 + 4t 2 (A-w) 2 



and A, 



_ 

4 1 + 4t 2 (A + co) 2 ' 
(«) _ 2r 2 (A + ^) 

" l + 4r 2 (A + o)) 2 ' 



+^( W ,i)+^^(a;,t) 



1 /A 



4 2t„ 



+ XA 



g-t/Tc 

-t/T C j. _ {2Ar c sin2(At + a) - cos2(At + a)} 



47or c 



- 1 



47or c 



(2Ar c sin 2a — cos 2a) 



1 

47or c 



(<") ' 



e -t/r c 

e * /Tc - {2Ar c sin2(Ai + a) - cos2(At + a)} 



47 T, 

(2Ar c sin 2a — cos 2a) 



4 3 W) = 



e t/Tc {sin 2 (At + a) + 2Ar c cos 2(At + a)} - {sin 2a + 2At c cos 2a} 



/ A") \ NO") 

J^\u,t) = + A4" ) -7T- (e~ t/2Tc (2t c (A - w) sin(A - uj)t - cos(A - uj)t) + l) 

V 2t c I 2t c I > 

+AA^ ) {e- t/2r -(sin[2a + (A - uj)t] + 2r c (A - uj) cos[2a + (A - u)t]) 



_ P ~t/2r c 



(sin(A — u>)t + 2r c (A — uj) cos(A — u))t) 



-(sin 2a + 2t c (A — uj) cos 2a) + 2t c (A — uj)} 



2t c 



{e _t/2Tc (2r c (A + w)sin[2a+ (A + u>)t] - cos[2a + (A + uj)t]) 



-(2r c (A + uj) sin 2a — cos 2a)}] 



4 5) (^o=(^+a4 w)N 



0") 



2t c 



{ e -t/2rc(2 Tc (A + w ) sin(A + w)t - cos(A + uj)t) + 1} 
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2t c 



.{ e -t/2rc(2 Tc (A - lo) sin[2a + (A - - cos[2a + (A - u)t]) 



— (2r c (A — w) sin 2a — cos 2a)} 

+A4 w) {e" t/2Tc (sin[2a + (A - w)t] + 2t c (A - w) cos[2a + (A - 

— (sin 2a + 2r c (A — w) cos 2a)} 

-AA^^e"*/ 2 ^ (sin(A + w)i + 2t c (A + w) cos(A + w)t) - 2t c (A + w)}] 



(A19) 



and 



4 6) ( W ,t)= a4 w) 



A 



2t c 



A 



2t> 



-{ e -* /2rc (sin(A - w)t + 2r c (A - w) cos(A - u)t) 



-2t c (A-w)} 

+ -^— {e-* /2Tc (sin[2a + (A + u)t] + 2t c (A + w) cos[2a + (A + u))i\) 
2r c 

— (sin 2a + 2r c (A + lo) cos 2a)} 

+\Af ) {e- t/2T "{2T c {\ + lo) sin[2a + (A + u)t] - cos[2a + (A + u))t]) 

— (2r c (A + w) sin 2a — cos 2a)} 

+\A^ ) {e- t ! 2T "{2 Tc {\ - w) sin(A - uj)t - cos(A - uj)t) + 1}] 



4 7 )( W) t) = [A4 w) A 



2t„ 



^{e- t/2Tc (sin[2a + (A - w)t] 



+2r c (A - w) cos[2a + (A-w)*]) - (sin 2a + 2r c (A - lo) cos 2a)} 
A {u) 

H — — {e"* /2Tc (sin(A + w)t + 2r c (A + w) cos(A + u)t) 
-2r c (A + w)} 

+AA^ } {e^* /2Tc (2r c (A - lo) sin[2a + (A - w)t] - cos[2a + (A - uo)t}) 
— (2t c (A — w) sin 2a — cos 2a)} 

+AA^ ) {e-*/ 2Tc (2r c (A + lo) sin(A + io)t - cos(A + w)t) + 1}] . 



(A20) 



(A21) 
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FIG. 1. Plot of o\ x {t) against time to show the short time behaviour of the variances for different temperatures with fixed 
parameters 70 = 1.0 and r c = 1.0. [ Inset : The same as in the main figure but for a higher temperature, ksT — 10.0 ] (units 
are arbitrary). 
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FIG. 2. Plot of cr xx (i) against time to show long time behaviour of the variances for different temperatures. Other 
parameters are same as in Fig.(l). [ Inset : The same as in the main figure but for a higher temperature, ksT = 10.0 ] (units 
are arbitrary). 

FIG. 3. Plot of a vv (t) against time to show long time behaviour of the variances for different temperatures. Other 
parameters are same as in Fig.(l) (units are arbitrary). 

FIG. 4. Plot of u xx (t) against time for different correlation times, r c with fixed parameters 70 = 1.0 and fcsT = 10.0 (units 
are arbitrary). 



FIG. 5. Same as in Fig. (4) but for fcsT = 0.0 (units are arbitrary). 

FIG. 6. Plot of u vv {t) against time to show long time behaviour due to vacuum fluctuations. Other parameters are same 
as in Fig.(l) (units are arbitrary). 

FIG. 7. Plot of quantum diffusion coefficient D q (t) against time for different temperatures and for 70 = 0.275 and r c = 1.0. 
[ Inset : Same as in the main figure but for a higher temperature fcsT = 10.0 ] (units are arbitrary). 
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